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Abstract 

We perform a complete Fourier analysis of the semi-discrete 1 — d wave equation obtained through a Pi dis- 
continuous Galerkin (DG) approximation of the continuous wave equation on an uniform grid. The resulting 
system exhibits the interaction of two types of components: a physical one and a spurious one, related to the 
possible discontinuities that the numerical solution allows. Each dispersion relation contains critical points where 
the corresponding group velocity vanishes. Following previous constructions, we rigorously build wave packets 
with arbitrarily small velocity of propagation concentrated either on the physical or on the spurious component. 
We also develop filtering mechanisms aimed at recovering the uniform velocity of propagation of the continuous 
solutions. Finally, some applications to numerical approximation issues of control problems are also presented. 

Resume 

On developpe une analyse de Fourier complete de l'equation des ondes unidimensionnelle semi-discretisee en 
espace obtenue dans 1' approximation numerique de l'equation des ondes par une methode de Galerkin discontinue 
(GD) Pi dans un maillage uniforme. On met en evidence la coexistence de deux composantes dans le systeme 
numerique : une physique, et une parasite liee aux discontinuites que la solution numerique permet. Chaque 
relation de dispersion contient des points critiques ou la vitesse de groupe correspondante s'annule. En suivant les 
constructions faites anterieurement pour le schema en differences finies, on construit d'une maniere rigoureuse des 
paquets d'ondes qui se propagent a une vitesse arbitrairement petite, concentres soit sur une composante ou sur 
l'autre. On developpe aussi des mecanismes de filtrage permettant de recuperer les proprietes de propagation des 
solutions de l'equation continue. Enfin, on presente une application a l'approximation numerique des problemes 
de controle. 
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Version frangaise abregee 

Dans cet article, on considere le probleme de Cauchy associe a l'equation des ondes unidimensionnelle 
semi-discretisee en espace obtenue dans l'approximation numerique de l'equation des ondes continue par 
une methode de Galerkin discontinue (GD) appelee Symmetric Interior Penalty Discontinuous Galcrkin 
(SIPG) (cf. [3]) dans un maillage uniforme de pas h en utilisant des polynomes du premier ordre. Ce 
schema produit deux relations de dispersion : une physique, associee a la partie continue de la solution 
numerique, ct une parasite, associee a ses sauts. On developpe d'abord l'analyse de Fourier de cette 
approximation qui met en evidence le fait que la vitesse de groupe correspondante a chaque relation de 
dispersion s'annule pour certains nombres d'onde. Cela nous permet d'adapter la construction de paquets 
d'ondes anterieurement realisee pour le schema en differences finics en [4], [9] et de construire des solutions 
numcriqucs concentrees sur l'un des deux modes possibles. Cela fait que la proprictc d'obscrvabilite bien 
connue pour l'equation des ondes continue, a savoir que pour un temps suffisamment grand l'energie totale 
des solutions peut etre estimee en fonction de l'energie localisee dans le complcmcntaire d'un ensemble 
compact (cf. [10]), n'est pas verifiee de maniere uniforme par rapport a h dans ce cas semi-discret. On 
exhibc ainsi un autre exemple de phenomene pathologique concernant les proprietes de propagation et 
dispersion des approximations numeriques classiques de l'equation d'ondes ([11], [4], [9]) et de Schrodinger 
(cf- [8]). 

Notre second objectif est de developper des mecanismes de filtrage pour construire des classes de 
donnecs initiales dans lesquelles la propriete d'obscrvabilite aie lieu uniformcmcnt par rapport a h. Pour 
ceci, on considere des donnees initiales avec des sauts nuls, assurant que l'energie totale des solutions soit 
dominee par l'energie des projections sur le mode physique, pour ensuite filtrer les hautes frequences par 
un algorithme bigrille ([6], [5], [7]). 

Nos resultats complement la litterature existante sur les methodes de Galerkin discontinus et en parti- 
culier [2], ou on demontre que la plupart de ces approximations sont spectralcment correctes, et [1], ou 
les proprietes dispersives et dissipatives de la version hp des semi-discretisations par GD de l'equation 
des ondes sont analysees. 

1. Fourier analysis of the Pi discontinuous Galerkin approximations of the 1 — d wave equa- 
tion. Based on discontinuous finite-element spaces, the discontinuous Galerkin (DG) methods can handle 
elements of various types and shapes, irregular non-matching grids and even varying polynomial order. A 
particular class of DG methods for elliptic and parabolic problems are the so-called interior penalty (IP) 
ones, where the continuity is weakly enforced across the element interfaces, by adding suitable bilinear 
forms, called numerical fluxes, to the classical variational formulations (see [3]). This note provides a 
further contribution to the analysis of DG methods. Among the existing and related works, we refer to 
[2], where the eigenproblcm associated to the Laplace operator discretized by means of DG methods is 
analyzed, showing that several DG methods provide spectrally correct approximations of the Laplace 
operator, and to [1], where the dispersive and dissipative properties of hp versions of various DG methods 
are studied. 

In this paper, we deal with the simplest setting of the Symmetric Interior Penalty Discontinuous 
Galerkin (SIPG) (cf. [3]) space semi-discretization of the 1 — d wave equation on an uniform grid Xj, j € 
Z, with first order polynomials. The numerical solutions given by the DG methods being discontinuous, 
their numerical representation in 1 — d consists in two values (possibly different ones) on every nodal 
point, representing the values to the left and to the right. An alternative, an often more convenient, 
way of representing the numerical solutions in the DG methods is in terms of the averages (denoted {•}) 
and the jumps (denoted [•]) along the interface, defined as {f}{x) = (f(x+) + f(x—))/2 and [f]{x) = 
f(x-) - f(x+). The finite element space is given by V h = U£ © Uj(, with U£ = span{0f, i e Z} n H^R) 
and U£ = span{0f , ieZ}, where 
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Figure 1. Basis functions for the Pi-DG methods: <f>^ (left) and </>^ (right). 
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Observe that (f>f are the typical basis functions of the Pi-classical finite element method, whereas </>/ 
are designed to represent the jumps at the nodal points. The functions {<p A , (fsf constitute a basis for 
V h . In this way, any / e V h can be uniquely represented as f(x) = Y^iez fi<t>t{ x ) + X«ez fi^i ( x ) = 
f A (x) + f J (x), where f A and f J are the continuous and the jump components of /, respectively. In this 
way, the piecewise linear discontinuous functions under consideration are perturbations of the classical 
piecewise linear and continuous ones by jumps added at each nodal point. 

Consider now the following Cauchy problem associated to the 1 — d continuous wave equation: 



d?u-dlu = 0, xeR, t>0, u{x,0) = u°(x), d t u(x, 0) = u 1 ^), x€ 



(1) 



which is well posed for (u ,^ 1 ) G ff^M) x L 2 (R), where H 1 ^) denotes the completion of C£°(M) with 
respect to the semi-norm || • ||^i (R) = \\d x ■ || L 2( K) . 

For each value s > 1 of the so-called penalty parameter, consider the symmetric bilinear form on Vh x Vh 

. . r r x i+i 

a s h (u,v) 



[/ 



u x (x)v x (x) dx - (M^Kt;,}^) + Mix^Mixj)) + -[u]( Xj )[v]( Xj ) 



and the following semi-discretization of the wave equation (1): 



u s h (x,t)eV h , df / u s h {x,t)v(x)dx + a s h (u s h (;t),v) = 0,VveV h . 



(2) 



The above variational problem (2) is completed by numerical initial data of the form u s h (x, 0) = (x) G 
Vh and d t u h (x,0) = u h (x) € Vh- The unknown u h (x,t), being an element of Vh for each t > 0, can be 

decomposed as u s h (x,t) = J^kez A k {t)4>£ ( x ) + Efcez Jk{t)^>l{x). Denoting by U h (t) = (A k (t), J k (t))' keZ , 
the problem (2) can be written in the following matrix form (' means transposition of a matrix): 

M h d*lf h {t) + R s h l} h {t) = Q, f h (0) = T} h '°, d t f h (0)=j} h -\ (3) 

where Mh and R s h are the block tri-diagonal symmetric mass and the stiffness matrices obtained by an 
infinite repetition of the stencils m/ l; r s hl where mh, r s h are the matrices 
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Set Hh = [— 7r/ft, 7r/ft]. For £ € 11^, let us denote by A h (£,i), ^(£,4) the semi-discrete Fourier trans- 
forms (SDFT) of the sequences of averages, A h (t), and of jumps, l h (t) (cf. [7]). Similarly, by A**®, 




Figure 2. Physical (black) and spurious (clotted black) dispersion relations, 1 (^) and X" p x (^), for s = 1.5 (top, left), 
s = 2 (top, right), s = 3 (bottom, left), s = 5 (bottom, right) compared to the ones corresponding to the continuous wave 
equation (blue), £, and to its finite difference (red) and Pi-classical finite element (green) semi-discretizations, and 
Ai(£). The marked points are wave numbers where the corresponding group velocities vanish. 

J^(£), i = 0,l, we denote the SDFTs of the initial data ~l h -\ ~1f h >\ Set U h (£,t) := (A h (^,t),J h (^t))'. 
The Fourier symbols of the mass and stiffness matrices are 

cos(gfe) ism(£h) \ /± sin 2(^> 
6 . . . m (ft = \ h 2 V 2 




2 - cos(gfe) ' - a - cos 2 (f) 

12 / V ° p 

Let S^CO := (M fc (f)) -1 .R£(f), (Af^))" 1 being the inverse of the matrix M h (£). The system (3) can 
be transformed into the following Cauchy problem associated to a system of two linear ODE's whose 
unknown is the vector function U (£,t), depending on the frequency parameter £: 

u&(z,t) + s' h (t)u h (t,t) = o, ^n hl *>o, u h (t,o) = u h >°(0, c/ t ' , (?,o) = ^ 1 (0^6% (4) 

Denote by A s ph<h (£), A s spih (£) the two eigenvalues of S£(f), by A^JO, h (0 their square roots and by 
Pph,h(Qi Ps P ,h(0 the two corresponding eigenvectors. The notation ph and sp stands for the "physical" 
and the "spurious" components, respectively. In this way, X* hh (t;), X s sp h (£) are the physical and the 
spurious dispersion relations. On the other hand, let Wft(£), Aft(£) be the dispersion relations corresponding 
to the finite difference (FD) and Pi classical finite element semi-discretizations of the 1 — d wave equation. 
The physical dispersion relation of the SIPG, X phh (£), satisfies: firstly, for h fixed, A* hf[ (^) — > Aft(£) as 
s —> oo, for all £ € Uh- On the other hand, for s > 1 and (Gil), fixed, X phh (£) —> |£| as h — > 0, which is 
the dispersion relation of the continuous wave equation. Contrarily, as s — > oo for fixed h or as h — > for 
fixed s > 1 and £ 6 lift, X s sp h (£) — > +oo. Moreover, for all s > 1, Ap h is an increasing function of £ 
and uJh((.) < Ap h < Aft(£), for all £ € lift, and the physical and spurious group velocities, d^X ph h (£) 
and d^X s ep h {£), have the following properties: 

1. For all s > 1, lim = 1 and hm 9 s A* p> ft(£) = 0. 

2. For all S G (1, oo") \ {3}, lim ^A^JO = t 1™ dtKp,h(0 = 0- 
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3- e lim a^ fc>h (0 = 1 and lim d £ X% >h (0 = -1. 

5— >±ir/h {— t±w/h 

4. For all s G (l,oo), the physical dispersion relation is increasing in £ for all s > 1, whereas the 
spurious one has several monotonicity ranges, according to the stabilization parameter s, as one can 
see in Figure 2. 

2. Localized waves and observability inequalities. For the continuous problem (1), it is well known 
that for any initial data (u ,^ 1 ) G x L 2 (M) and any observability time T > 2, there exists a 

constant C(T) > s.t. the following observability inequality holds (cf. [10]): 

T 

E{u\u x ) sC C(T) J E n (u°,u\t)dt, (5) 
o 

where Q := R \ (—1,1) and E(u®,v}) and Eq(u° 1 u 1 ,t) denote respectively the total energy (which is 
conserved in time) and the energy concentrated in Q at time t, given explicitly by 

Eiu^u 1 ) = ^ J {\d t u{x,t)\ 2 + \d x u{x,t)\ 2 )dx, E n (u°,u\t) = ^J(\d t u(x,t)\ 2 + \d x u(x,t)\ 2 )dx. 
r n 

The time T* = 2 is sharp, given by the so-called Geometric Control Condition (GCC), requiring all 
rays of Geometric Optics to enter the observation region during the observability time. When the GCC 
does not hold, the observability property fails because of the Gaussian beam solutions localized around a 
bi-characteristic ray escaping the observation region (see [4]). 

We also analyze the SIPG version of the observability inequality: 

T 

E'tfhflrfh,!) < C s {T) J E s nth {f h >\f h >\t)dt. (6) 



Here, E s h (ff h '°, C^ 1 ) is the discrete total energy, conserved in time, given by 
E s h {f h '°,lf h ' 1 ) = ^(R s jf h '°,f h '°) + {M h tf h ' 1 ,t h ' 1 \ 

where (•, •) is the inner product in £ 2 , and E^ h (u h '°,u h ' 1 ,t) is the discrete energy concentrated in il 
at time t, defined as above but with (•, •) replaced by the local scalar product (•, which is the inner 
product in ? 2 ({j : Xj G SI}). For all T, h > 0, the inequality (6) holds, with a finite constant C^(T). Our 
goal here is to analyze its behavior as h — > 0. 

When the vector valued initial data U h ' 1 in (4), i = 0, 1, are of the form 

u h H0 - PUAO^iO, (7) 

the corresponding solutions of (4) involve only the physical dispersion relation: 

Considering solutions concentrated in wave packets and a stationary phase like argument allow showing 
that, whatever s and T are, C^(T) blows up at an arbitrarily large polynomial rate as h — > 0: 
Proposition 0.1 Let T > be given with a semi-discrete ray x ph (t) = x* — td^X ph x (?7o) that does not 
enter the observation region in time T. Consider 7 := 7(/i) > such that 7 >> 1 and h'y << 1. For 
<t> G 5(R), consider the semi-discrete wave equation (4) with initial data U h,l {^) satisfying (7) with 

u h '°(0 = (^) cxpHx*(£ - eo))Xn,(0 and g^tf) = iX ph , h (0^°(0- (9) 



Then for all a e R + , the observability constant C£(T) in (6) satisfies C£(T) ^ C a (4>,T,s)^ a . 
A similar result holds when the numerical solution is concentrated on the spurious mode. 

3. Filtering mechanisms. In what follows, we introduce a filtering mechanism aimed at recovering the 
uniformity as h tends to zero of the observability constant C^(T) in (6) within a subclass of solutions of 
the numerical approximation scheme. For 5 e (0, 1), set 11^ := [—Tr5/h,ir5/h} and let us define the space 

of Fourier filtered data J| = { / € £ 2 (hZ) : supp(/ ?l ) C It can be proved that if in (3) we consider 
initial data lf h ^, i = 0,1, verifying the condition (7) and such that u h >' 1 e 1^ for i = 0,1, then there 
exists an uniform time T*£ such that, for all T > T*£ , the observability inequality (6) holds uniformly as 
h — >• 0. These data lead to solutions whose energy is concentrated on the low frequencies of the physical 
mode for which the group velocity of propagation is uniform. 

However, our goal (as described in the pioneering work [5]) is to introduce a filtering mechanism that 
does not require the use of the Fourier transform, but that rather might be implemented in the numerical 
mesh, directly. The filtering mechanism we propose can be implemented in two steps. First, the initial 
data are taken so that their jumps vanish, i.e. ~f h ^ = 0,i = 0, 1. Second, their average part, ~A htt , is 
obtained by a bi-grid algorithm (analyzed in [6]), i.e. A\j = (^sy+i + ^2j-i)/2, Vj e Z, \/i = 0, 1. 

Although the corresponding solutions also excite the spurious spectral component, their energy is 
concentrated on the low frequency physical components. Thus, using the arguments in [7], one can show 
that, for T large enough independent of h, the semi-discrete observability inequality (6) holds uniformly 
in this class of filtered numerical solutions too. By duality, this implies a result of uniform (with respect 
to h) controllability of a suitable projection of the solutions of the numerical approximation scheme. 
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